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Abstract 

Protocols of quantum energy teleportation (QET), while retaining 
causality and local energy conservation, enable the transportation of 
energy from a subsystem of a many-body quantum system to a distant 
subsystem by local operations and classical communication through 
ground-state entanglement. We prove two energy-entanglement in- 
equalities for a minimal QET model. These relations help us to gain a 
profound understanding of entanglement itself as a physical resource 
by relating entanglement to energy as an evident physical resource. 



1 Introduction 



The relationship between energy and information has been investigated 
extensively in the context of computation energy cost including a modern 
analysis of Maxwell's demon pP-[2]- In this Letter, we show a new energy- 
information relation from a different point of view. Recently, it has been 
reported that energy can be transported by local operations and classical 
communication while retaining local energy conservation and without break- 
ing causality |3]-|5]. Such protocols are called quantum energy teleportation 
(QET) and are based on ground-state entanglement of many-body quantum 
systems including spin chains [3J, cold trapped ions [3] and quantum fields [3]. 
By performing a local measurement on a subsystem A of a many-body sys- 
tem in the ground state, information about the quantum fluctuation of A can 
be extracted. Because the post-measurement state is not the ground state in 
general, some amount of energy is infused into A as QET energy input during 
this measurement, and the ground-state entanglement gets partially broken. 
Next, the measurement result is announced to another subsystem B of the 
many-body system at a speed much faster than the diffusion velocity of the 
energy infused by the measurement. Soon after the information arrives at B, 
energy can be extracted from B as QET energy output by performing a local 
operation on B dependent on the announced measurement data. The root 
of the protocols is a correlation between the measurement information of A 
and the quantum fluctuation of B via the ground-state entanglement. Due to 
the correlation, we are able to estimate the quantum fluctuation of B based 
on the announced information from A and devise a strategy to control the 
fluctuation of B. By the above-mentioned selected local operation on B, the 
fluctuation of B can be more suppressed than that of the ground state, yield- 
ing negative energy density around B in the many-body system. The concept 
of negative energy density has been investigated in quantum field theory for a 
long time [6]. Quantum interference among total energy eigenstates can pro- 
duce various states containing regions of negative energy density, although 
the total energy remains nonnegative. The regions of negative energy density 
can appear in general many-body quantum systems by fixing the origin of 
the energy density such that the expectational value vanishes for the ground 
state. In spite of the emergence of negative energy density, the total energy 
also remains nonnegative for the general cases. In the QET protocols, dur- 
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ing the generation of negative energy density at B, surplus positive energy 
is transferred from B to external systems and can be harnessed as the QET 
output energy. Here it should be emphasized that this output energy existed 
not at A but at B even before the start of the protocol and was hidden in- 
side the zero-point fluctuation of B. Of course, this zero-point energy is not 
available by usual local operations for B. However, by using a local operation 
dependent on A's information, it becomes possible to dig out B's zero-point 
energy by pair creation of the positive output energy from B and the negative 
energy of B. Hence, we do not need to hire any physical carrier of energy from 
A to B like electric currents and photons, at least, during short-time QET 
processes. Needless to say, after the completion of QET process, the positive 
energy of A compensates for the negative energy of B during late-time free 
evolution of the many-body system. The amount of output energy from B is 
upper bounded by the amount of input energy to A. 

Using the usual protocols of quantum teleportation, quantum states can 
be transported from one party to another by the consumption of shared 
entanglement between the two parties [7]. As is well known [S], transfer 
of a large number of quantum states requires a large amount of consump- 
tion of shared entanglement as a physical resource. Taking into account the 
fact, it seems natural for the QET protocols to expect that a large amount 
of teleported energy also requests a large amount of consumption of the 
ground-state entanglement between A and B. If such a non-trivial relation 
exists between teleported energy and breaking of ground-state entanglement 
by measurement, the relation may shed new light on the interplay between 
quantum physics and quantum information theory. In this Letter, the first 
example of the energy-entanglement relation for a minimal QET model is pre- 
sented. The minimal QET model is the smallest physical system for which 
non-trivial QET can be implemented; this model consists of two qubits with 
an interaction of the Ising spin chain in the presence of a transverse mag- 
netic field. We explicitly show that for the minimal model, the consumption 
of entanglement between A and B during the measurement of A is lower 
bounded by a positive value that is proportional to the maximum amount 
of energy teleported from A to B. In addition, we obtain another inequality 
in which the maximum amount of energy teleported from A to B is lower 
bounded by a different positive value that is proportional to the amount of 
entanglement breaking between A and B by the measurement of A. These 
energy-entanglement inequalities are of importance because they help in gain- 
ing a profound understanding of entanglement itself as a physical resource 
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by relating entanglement to energy as an evident physical resource. 



2 Minimal QET Model 



First of all, we introduce the minimal QET model. The system consists of 
two qubits A and B. Its Hamiltonian is the same as that of the Ising spin chain 
in the presence of a transverse magnetic field as follows: H = Ha + Hb + V, 
where each contribution is given by 

and h and k are positive constants with energy dimensions, (j\ (cr^) is the 
x-component of the Pauh operators for the qubit A (B), and a\ (cr^) is the 
z-component for the qubit A (B). The constant terms in Eqs. (II])-® are 
added in order to make the expectational value of each operator zero for the 
ground state \g): 

{g\HA\g) = {g\HB\g) = {g\V\g) = Q. 

Because the lowest eigenvalue of the total Hamiltonian H is zero, if is a 
nonnegative operator: H > 0. Meanwhile, it should be noticed that Hb has 
a negative eigenvalue, which can yield negative energy density. The ground 
state is given by 





where |±)a is the eigenstate of (o"b) with eigenvalue ±1. 

When local operations are performed to this system, we connect the sys- 
tem with a local appratus and (or) add local external forces like magnetic 
field. However, time interval of the operations is assumed to be quite short. 
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as often argued in quantum information theory. Thus interaction terms be- 
tween the system and the external apparatus vanish in the Hamiltonian after 
the short-time operation. 

Let Sma denote a set of POVM measurements [8] for A which measure- 
ment operators Ma{^j) with measurement output /x commute with the in- 
teraction Hamiltonian V . The measurement operator M^(yu) takes the form 
of 

MA(/i) = e^'- K + e*"-/^cr^) . 
The coefficients m^, and 5^ are real constants which satisfy 

^m^/^ cosa^ = 0. 
The POVM corresponding to Ma{ij) is defined by 

nA(/i) = MA(/i)^MA(/i), 

which satisfies the completeness relation, 

^n^(/i) = iA. 

By introducing the emergence probability Pa{^J') of output /i for the ground 
state and a real parameter qAifJ-), the POVM is written as follows: 

By taking suitable values of m^, l^, and a^, all values of PA^fJ') and gyilAt) 
are permissible as long as they satisfy J2fj.PA{fi') = 1, = and 

> \(lA{fJ')\- The post-measurement state of the two qubits with output 
/i is given by 

\Ai^,)) = ^l=MAi^^)\9). 

This measurement excites the system and the average post-measurement 
state has a positive expectational value Ea of H, which energy distribution 
is localized at A. In fact, the value defined by 

EA = J2{g\MAi^yHMA{fx)\g) 
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is computed as 



This infused energy Ea is regarded as the QET energy input via the measure- 
ment of A. During the measurement, Ea is transfered from external systems 
including the measurement device with a battery respecting local energy con- 
servation. The QET energy conservation law during local measurements to 
a qubit of a spin chain has been discussed in [3]. 

The key feature of this model is that any measurement of Sma does not 
increase the average energy of B at all. By explicit calculations, the average 
values of the Hamiltonian contributions Hb and V are found to remain zero 
after the measurement and are the same as those of the ground state. This 
measurement does not yield instantaneous change of V. Therefore we have 
no direct force of A affecting B after the measurement. Thus, we cannot 
extract energy from B by the standard way soon after the measurement. 
In fact, if any local unitary operation Wb independent of A's measurement 
result is performed to B, the post-operation state u is given by 



Wb \^MA{fi)\g){g\MA{fiyj Wl 

The energy difference after the operation is calculated as 

Tr [uH] -Ea = {g\Wl {Hb + V) WB\g), (5) 
where we have used 

W^HaWb = HaW^Wb = Ha, 

wI{Hb + V)Wb, M^(/i)] =0, 
and the completeness relation of the POVM's: 

J2Ma{iiVMa{i^) = 1a. 

From Eq. (|5]), it is proven that the energy difference is nonnegative: 
Tr [uH] -Ea= {g\WlHWB\g) > 0, 
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because of a relation such that {g\W^HAWB\g) = {g\HA\g) = and the 
nonnegativity of H. Therefore, as a natural result, no local operation to B 
independent of /i extracts energy from the system. 

After a while, the infused energy Ea diffuses to B. The time evolution of 
the expectational values Hb and V of the average post-measurement state is 
calculated as 

and (y{t)) = 0. Therefore, energy can be extracted from B after a diffusion 
time scale of this is just a usual energy transportation from A to B. 
Amazingly, the QET protocol can transport energy from A to B in a time 
scale much shorter than that of the usual transportation. In the protocol, 
the measurement output fi is announced to B. Because the model is non- 
relativistic, the propagation speed of the announced output can be much 
faster than the diffusion speed of the infused energy and can be approximated 
as infinity. Soon after the arrival of the output fi, we perform a local operation 
Ub{ij) on B dependent on /x. Then, the average state after the operation is 
given by 

p = 5^f/B(/i)M^(/i)|^7)(^7|M^(/i)tf/^(/i)t. 

In Figure 1, a schematic diagram of this QET model is presented. The 
expectational value of the total energy after the operation is given by 

Tr[pi7] = Y,'^9\MA{^^)^UB{^^)^HUB{^^)MA{^^)\g). 

On the basis of the fact that Ub{ij) commutes with Ha and Eq. Eb is 
computed as 

Eb = Ea- Tr [pH] = -Tr [p{Hb + V)] . 

Further, on the basis of the fact that M^(/i) commutes with UB{fi), Hb and 
V, the energy can be written as 

Eb = - $^(^|nA(/x) (HBi^i) + V{^)) \g), 
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where the /i-dependent operators are given by Hsifi) = Ub{ij)'^ HbUb{ij) 
and V{ij) = Ub{ij)'^VUb{ij)- Here, let us write the general form of Ub{ij) as 
follows: 

Ub^Ii) = cosu^ + in^ ■ afisinc;^, 

where uj^ is a real parameter, = {rix^, Uy^, riz^) is a three-dimensional unit 
real vector and is the Pauli spin vector operator of B. Then, an explicit 
evaluation of Eb becomes possible. The result is expressed as 

where is given by 

Q{lj) = cos (2a;^) - hkqA{fi)nyf, sin(2a;^) - X{fi), (6) 

where X{fi) is defined by 

^(/^) = Pa(/^) [h^ (1 - nl^) + 2A;2 (l - nl^)] - 3hkqA{fi)n^f,n,f,. 

In order to maximize the teleported energy Eb for a given POVM measure- 
ment of A, let us first maximize Q{fi) in Eq. (j6]) by changing the parameter 
Ufj^. This maximum value is calculated as 

maxQ(/i) = \/x{fiy + [hkqA{fi)ny^f - X{n). (7) 

Next, let us introduce a parametrization of n^^ and n^^ as 71^^ = -^/i cos 
and = -y/i sin for fixed z = l—riy^ which runs over [0, 1], where ip^ is a 
real parameter. It is observed that max^^ in Eq- (E]) is a monotonically 
decreasing function of X{^). Thus, we must find the minimum value of X{n) 
in terms of ip^. By using the parametrization, we can minimize X{^) as 

rninX(/x) = (l - {h' + 2^)-^ {h^ - 2P)'p^(/x)2 + gh^k^qAi^i)^. 

Therefore, the maximum value of max^^^ in terms of ipf^ is obtained as 
follows: 



max = \ I minX(/i) I + K^k'^qA^nYil — z) — minX(/i). 



Next, in order to maximize max^^^^^^ Q(a*) in terms of let us write it as a 
function T{z) of z: 



T(z) = m.axQ{jj,) = y (a — bz) + c{l — z) — {a — hz) , 
where a, h and c are positive constants given by 

b = ^ {h' + 2^) + y {h? - 2k^fpAW + ^h?k^A{li)\ 

The derivative of T{z) can be calculated as 

2^{a-hzf + c{l- z) 
where t{z) is a function given by 

t{z) = -c + 26 {a-hzf + c{l- z) - (a - 6^)^ . 

It can be verified that t{z) and dzT{z) are nonpositive for z G [0,1]. This 
verification can be done as follows. Let us first consider an equation t{z) = 0. 
It turns out that, in the transformation of this equation for solving z, the 
dependence of z gets lost and we get just a constraint condition on Pa(a*) and 
gA(/i) such that qA^fJ')^ {PAif^Y - qA{lA'^) = 0. Thus, if pAifJ'Y = qA{^)'^ or 
qAifJ')'^ = 0, the equation t{z) = holds for all z G [0, 1]. If pAi^^f 7^ qA{f^)'^ 
and qAifJ'Y 7^ 0, the solution z does not exist and t{z) has a definite sign 
for z G [0, 1]. In order to check the sign, let us substitute z = 1 into t{z). 
Then, when a > 6, we get t (1) — —h?k'^qA{lJif' < 0, and when a < 6, t (1) = 
—'&h?k'^{pA{lj)'^ — qAif^Y) < 0. Thus, it is verified that t{z) and dzT{z) are 
nonpositive. Therefore, T{z) takes the maximum value at 2; = 0. This implies 
that Q{fJ^) can be maximized as maxt/g(^) := max^^^^^^^ (^(/u) = T(0). 

This leads to our final expression of the maximum teleported energy for the 
measurement, which is clearly nonnegative, as follows: 



max Eb 

Ub{iA 



+ 2k^ 



^PA{li) 



{h^ + 2k^f PAii^y 




(8) 
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The operation C/max(A*) which attains the maximum of teleported energy is 
given by 

t^max(/^) = COS Qf, + lu'^ sin f]^, 
where is a real constant which satisfies 



cos (2Qju) 
sin (20^) 



hkqAjjj) 



Besides, the teleported energy can be maximized among POVM measure- 
ments of Sma- This is achieved when each POVMs are proportional to pro- 
jective operators and given by 



+ 2k^ 

max Eb = , - 



'1 + 



{h? + 2F)^ 



- 1 



3 Relation between Entanglement Breaking 
and Teleported Energy 

Next, we analyze entanglement breaking by the POVM measurement of A 
and show two inequalities between the maximum teleported energy and the 
entanglement breaking. We adopt entropy of entanglement as a quantitative 
measure of entanglement. The entropy of a pure state I^'as) of A and B is 
defined as 

Tr[|*^B)(*^B|]lnTr[|*^B)(*^B|] ■ 

Before the measurement, the total system is prepared to be in the ground 
state \g). The reduced state of B is given by 

PB = Tj[\g){g\]. 



Sab^-"^ 
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After the POVM measurement outputting the state is transferred into a 
pure state The reduced post- measurement state of B is calculated as 

Pb{i^) = ^^Tr[UA{i^)\g){g\]. 

The entropy of entanglement of the ground state is given by 

Sab{9) = - Tr \pb Inps] 

and that of the post-measurement state with output // is given by 

SAsii^) = - Tr lnps(/i)] . 

By using these results, we define the consumption of ground-state entangle- 
ment by the measurement as the difference between the ground-state entan- 
glement and the averaged post- measurement-state entanglement: 

H^Sab = Sab{9) - ^Pa{ij)Sab{ij)- 

Interestingly, this quantity is tied to the quantum mutual information be- 
tween the measurement result of A and the post-mcasTircmcnt state of B. 
Let us introduce a Hilbert space for a measurement pointer system A of the 
POVM measurement, which is spanned by orthonormal states corre- 
sponding to the output 11 satisfying = Then, the average state 
of A and B after the measurement is given by 

^AB = ^pa{ij)\iJ'a){iJ'a\ ®pb{ij)- 

By using the reduced operators $^ = Tr_B [$^5] and $b = Tr^[$^5], the 
mutual information Iab is defined as 

Iab^-"^ [^A In - ^ [^B In + IV [^ab In ^ab] ■ 

A B AB 

By using Tr^ [^ab] = Y.,,Pa{p)\pa){pa\ and Tr^ [^ab] = T.hPa{h)pb{h) = 
Pb, it can be straightforwardly proven that Iab is equal to ASab- This 
relation provides another physical interpretation of ASab- 

Next, let us calculate ASab exphcitly. All the eigenvalues of Pb{p) are 
given by 
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1 ± 1 / cos^ <j + sin <^ 



(9) 



where <^ is a real constant which satisfies 

h 



cos<^ 



sin<^ 



k 



The eigenvalues of ps are obtained by substituting g^l/^) = into Eq. ([9]) 
By using Xs{p), ^Sab can be evaluated as 



^Sab = ^PA{fJ')fl 



PA{pf 



(10) 



where fi{x) is a monotonically increasing function of x G [0,1] and is defined 

by 



//(x) = ^ (l + V cos^ + x sin^ j In Q ^1 
- ^1 — \/ cos^ + X sin^ ^ j In ^- ^1 



+ 



(1 + cosg) In I - (1 + cos?" 



cos^ ? + X sin q 



cos^ ? + X sin ? 



(1 — cos?") In I - (1 — cos? 



It is worth noting that the optimal teleported energy max[/^(^) Eb in Eq. ([8]) 
takes a form similar to Eq. (fTOjl as 



max Es = ^^^(/i)/^ ( 
t^s(w — V 



PA(/i)V ' 



11) 



where fE{x) is a monotonically increasing function of x G [0, 1] and is defined 
by 



/s(x) = VWTk^ {1 + sm\) 



'1 + 



cos^ ?" sin^ ?" 



(l + sin'?)' 

Expanding both //(x) and /£;(x) around x = yields 



X — 1 
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fi{x) = In X + O(x'), 

4 cos ? 1 — cos ^ 

X / N / , 9 , , n COS^^Sin^? , ^, 2^ 

2 (1 + sm 

By deleting x in the above two equations, we obtain the following relation 
for weak measurements with infinitesimally small qA{f^)'- 

. l + sin\ 1 + cos^max^^(^)EB ^ 

2 cos^ ^ 1 - cos <^ V/;,2 + 

It is of great significance that this relation can be extended as the following 
inequality for general measurements of Sma '■ 

> In l + cos.max^,(,)i?^ 

- 2cos3^ 1-cos^ yFTP 
This inequality is one of the main results of this Letter and implies that a large 
amount of teleported energy really requests a large amount of consumption 
of the ground-state entanglement between A and B. The proof of Eq. ( IT2ll 
is as follows. Let us introduce two rescaled functions as follows. 

J7 , . COS^ / 1 + cos^\ , , 



sin g \ 1 — cos 

1 2(l + sin\) ^^^ 

fE{x) = ^ — r^fEix). 

V + cos^ q sm q 



It can be easily shown that /^'(x) is a convex function: 

dlfEix) < 0. 

From this convexity and = x + O(x^), the function fE{x) satisfies 

hix) < X (13) 

for < X < 1. On the other hand, it is observed that fi{x) is a concave 
function, as shown below. The derivative of fi{x) is computed as 



1— cos f 
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where y{x) = a/ cos^ <; + x sin^ q and gi{y) is a positive function of y defined 
as ^ ^ 

9i{y) = - In- -. 

y 1 

It should be noted that y{x) is a monotonically increasing function of x. The 
derivative of gi{y) is calculated as 



dygiiy) = 

where si{y) is a function of ?/ given by 
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y 



si{y) = ^ - In 

and satisfies a boundary condition as s/(0) = 0. It is also easy to show that 
the derivative of si{y) is positive for y > 0: 

dysiiy) = 7-2 > 0. 

(1 -r) 

Thus, sj{y) and dygj{y) are positive for y G [cos?", 1]. From these results, it 
has been proven that //(x) is concave for x G [0, 1]: d1fj{x) > 0. From this 
concavity and fi{x) = x + 0{x^), it is shown that //(x) satisfies the relation 

fi{x) > X. (14) 

Because of Eqs. ffT51) and ffTll) . we obtain the following inequality: 

fi{x) >x> Je^x). 

This implies that the energy-entanglement inequality in Eq. (|T2|) holds. In 
addition, we can prove another inequality between energy and entanglement 
breaking. Because the convex function /^(x) and the concave function //(x) 
are monotonically increasing functions of x G [0, 1] which satisfy /^lO) = 
//(O) = 0, we have the following relation: 

fE{x) > ^fl{x). 
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Consequently, the following inequality, which is another main result of this 
Letter, is obtained for all measurements of Sma'- 



^ ^ + fc2 Wa - 3 cos2 ^ - 2 + cos2 d 

max_fcB>- — — — r — — 7 — ^ — ^iXbAB- (15) 

- (l + cos,)ln(^) + (l-cos,)ln(^) 

This ensures that if we have consumption of ground-state entanglement 
IS.Sab for a measurement of Smai can in principle teleport energy from A 
to B, where the energy amount is greater than the value of the right-hand- 
side term of Eq. ( !T5|) . This bound is achieved for non-zero energy transfer by 
measurements with = The inequalities in Eq. ( !T2|) and Eq. 

(fTSjl help us to gain a deeper understanding of entanglement as a physical 
resource because they show that the entanglement decrease by the measure- 
ment of A is directly related to the increase of the available energy at B as 
an evident physical resource. 
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Figure Caption 



Figure 1: Schematic diagram of the minimal QET model. A POVM mea- 
surement is performed on A with infusion of energy Ea- The measurement 
result fjL is announced to B through a classical channel. After the arrival of 
II, a unitary operation dependent on /i is performed on B with extraction of 
energy Eb- 
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